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ABSTRACT 


The effects of horizontal refractivity gradients on the accuracy of 
laser ranging systems were investigated by ray-tracing through three- 
dimensional refractivity profiles. The profiles were generated by 
performing a multiple regression on measurements from seven or eight 
radiosondes, using a refractivity model which provided for both linear and 
quadratic variations in the horizontal direction. The range correction 
due to horizontal gradients was found to be an approximately sinusoidal 
function of azimuth having a minimum near 0“ azimuth and a maximum near 
180® azimuth. The peak-to-peak variation was approximately 5 centimeters 
at 10“ elevation and decreased to less than 1 millimeter at 80“ elevation. 
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1. INTRODUCTION 


1 


Laser ranging can be used to precisely measure the distance from a 
point on earth to an orbiting satellite. Ranging accuracy is limited by 
atmospheric refraction and scattering. In this paper, the effects of 
atmospheric refraction will be Investigated. Atmospheric refraction will 
Increase the optical path length to an orbiting satellite by over 13 meters 
when the satellite Is at 10° elevation, and by about 2 meters when the satel- 
lite Is directly above the tracking station. 

Numerous formulas have been developed to estimate the range error due to 
refraction. These formulas use surface measurements of pressure, temperature, 
and relative humidity to predict the range error. The analysis developed 
by Marini and Murray [1] is used in this paper to derive a new surface 
correction formula. The formula is derived under the assumption that the 
atmospheric refractivity is spherically symmetric, that is, the index of 
refraction is a function of height only. Horizontal refractivity gradients 
will introduce errors into the surface correction formula. Zanter, Gardner 
and Rao[2] investigated this effect by ray tracing through a refractivity 
model having a linear dependence on horizontal position. In this report, 
a quadratic refractivity model is used to determine the accuracy of the 
range correction formula by ray tracing at various azimuth and elevation 
angles . 









2. SLRFACE CORRECTION FORMULA 


The geometry of the laser ranging problem is shown in Fig. 1. The 
laser pulse travels a curved path from the tracking station (H meters above 
sea level and r^ meters from the center of tiie earth) to the satellite (r^^ 
meters from the center of the earth). is the arrival angle of the 

pulse, while E is the satellite's actual elevation angle. 6 is the angle 
of the ray with the horizontal (a function of height). The problem is to 
find R, the straight-line distance to the satellite by measuring the curved 

distance R . 

o 

The derivation of a surface correction formula requires the index of 
refraction n and the corresponding phase refractivity N [1]: 


10^ (n 


1 ) 


= 287.604 + 


1.6288 0.0136 1 i 273.15' ' P 

^2 ^4 111013.25/ T - 
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273.15 

U013.25; 


where X = wavelength of laser in ml -rens 

P = atmospheric pressure in mill J oars 

T = temperature in degrees Kelvin 

e = partial pressure of water vapor in millibars. 


The group refractivity N and group index of refraction n must also be 

8 8 

known [ 1 ] 
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where f(X) = 0.9650 + n.'')164/>^ 


0.000228/ 
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( 2 - 1 ) 


According to the theory of geometnh 


iptii.s, the optical path Length 





n dl 
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( 2 - 2 ) 


where dl Is an incremental length along the ray path C. If the group Index 
of refraction is a function of height only (spherical symmetry), the ray 
path will lie entirely in a plane, and 

dr 


dl = 


sin (0) 


so (2-2) becomes 
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sin (6) 


dr 
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r’^l N 

J. sin 
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dr + 




dr 


sin (6) 


^0 


The difference between the curved and straight-line distances is the 
range error AB 

r~ 

r. 


AR = R - R = 10 
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-6 
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A. 


sin (6) 


dr + 
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f ^ - 4£ 

j So.1- (0) 

. 0 


- R 


(2-3) 


The bracketed term corresponds to the geometrical error in path length 
and has been evaluated in [!]• 
r. 


dr 


sin (0) 


- R == 
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sin^ (0 q) 


i 10-^2 


dh 


(2-4) 
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The first term of (2-3) corresponds to the velocity error. Since the satel- 
lite is above the atmosphere (r^ > 70 km) where = 0, ttie upper limit can 


be extended to infinity and a change of variable from r to h gives 

r. 
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sin (0) 
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sin (9) 


dh 


(2-5) 
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To evaluate the right-^hand side of (2-5), Marini and Murray expanded 
the Integrand in inverse powers of sin (6 q)* 0q can be related to 0 by 
Snell's law for spherically symmetric media: 

nr cos (0) » * 

This formula is solved for - v 

sin (0) 


sin (0) sin (0^) 


i.fVo) 

1 + 


TTn * 


sin (0) 


tan (0 q) 


can be expanded using approximations detailed in Appendix A, 


sin (0) ” sin (0 q) 


^ - 10’^ (Nq - n) I 

1 - 5 + “ 

tan (0 q) 


^ - 10"® (N - N) 
0 
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tan 0, 


'(2-6) 


Using (2-5) and (2-6) in (2-3), we obtain 
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(2-7) 


Marini and Murray express ;,.ieir result in terms of the actual satellite 
eleva*’ion angle E rather than t! arrival angle 6^. The conversion can be 


Unless otherwise noted, the limits on all integrals in this report will 
be from 0 to 
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made using tSe first term of the angular correction, il]: 

e - E = lO'^N- cot (E) . 

0 0 


Adding E to both sides and taking the sine 

(e + 10~^ 


sin (0 q) = sin 
..- 6 . 


Nq cot (E)j 


Nq “ 300, so for E ^ 10°, 10 Nq cot (E) is small, and 

sin (0q) = sin (E) + 10 ^Nq cot (E) cos (E) 
^ ^ 


sin (6q) - sin (E) ^^^2 ^ 


Now 10 ^Nq cos^ (E) << sin (E) for E ^ 10°, giving 


^ 10"Sq[ 1 - sin^ (E)] 


sin (9q) sin (E) 


sin^ (E) 


10"^N, 


sin (E) . 3 . . * 

sin (E) 

Applying (2-8) to (2-7), we obtain 
sin (E) 
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Many of the terms in (2-9) can be igrored, as they either make a con- 
tribution of less than 1 mm to AR at E * 10°, or they cancel each other. 
Making these simplifications, which are detailed in Appendix B, 
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The underlined terms in (2-10) appear in the correction formula developed 
by Marini and Murray [1] and have been evaluated by them. The last integral 
in (2-10) was not used by Marini and Murray and is evaluated in Appendix C. 

The integral evaluations, taken from Marini and Murray [1] and Appendix C. are 
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2 2 
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where f(X) * 0.9650 + = — + ; 

,2 ,4 

A X 


F(e,H) = 1 + 0.0026 COL (20) - 0.00031H 


( 2 - 11 ) 


K = 63 + 0.00968 cos (20) - 0.001041^ + 0.00G'’435P 
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Latitude of laser site 



H “ surface height at laser site (in km) 


H 


i 



P - surface pressure at laser site (in mb) 
s 

■ surface te^>erature at laser site (in ®K) 

e - surface partial pressure of water vapor at laser site (in irf)) 
s 

\ = wavelength of laser (in microns). 

Equation (2-10) can be expressed in a continued fraction form, similar to 
that u^ed by Marini and Murray [1] 


AB = f(X) 


sin (E) 


3 5 

Jin (E) sin (E) 


= f.A) 


sin (E) + 


sin (EX + 


B/A 

C/B 

sin (E) + 0.17 


( 2 - 12 ) 


where A 


I0.002357P + O.OOOlAle ] + 1.0842 x 10~®P T K 
Fv8»H/ S s s s 

p2 

- 9.4682 X 10~® ~ 

" p2 

B » 1.0842 X 10"®P T K + 4.7343 x lo'® ^ i - \ 

s s T 3 - 1/K 

s 

2 

-1 7 K 

C = 1.4961 X 10 ^^p T 


s s 2 - K 


The 0.17 is an empirical constant which compensates for approximations 
made in the derivation of (2-12). The surface correction formula developed 
by Marini and Murray [1] contains an empirical constant (0.01) which 

c/b 

replaces the tern gin (F.) + 0 17 ^2-12) (see Appendix D) . 





Marini and Ibirray compared their fomila with a range correction 
obtained by ray tracing through a spherically syasietric ataosphere. The 
assuaqption o£ spherical synnetry means that the refractivity is a functicm 
of height only (independent of horizontal position). Any nonsyosetric 
behavior of the ataosphere would introduce errors in a ray trace made under 
> the spherical syaaetry assumption. 

To investigate the errors introduced by this assus^tion, Zanter, 
Gardner and Rao [2] assumed that the refractivity at a given height had 
a linear dependence on position: 

N - N + 0N- + ♦ sin (6) N. (3-1) 

r (7 9 

where M ■ refractivity 

9 -> colatitude ~ 90* - latitude 
9 ■ longitude 

N ■ coefficients to be determined for each height. 

r u 9 

6 is proportional to horizontal displaceaent in the north-south direction, 
while 9 sin 6 is proportional to horizontal displacement in the east-mst 
direction. The three coefficients can be determined at any particular 
altitude from a knowledge of the refractivity at a minimum of three 
points at that altitude. 

Refractivity does not exactly follow the linear model of (3-1), and 
a more accurate quadratic model was used for this investigation: 

N - N + 0N„ + 9 sin (0) + 09 sin (0) + 0^N„. + 9^ sin^ (9) . 

r O 9 09 Ho 9 r 

(3-2) 

The six coefficients can be determined at a particular altitude from a 
minimum of six refractivity measurements. 
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Refractlvity is calculated from radiosonde measurements of pressure, 
temperature, and relative humidity. Errors in radiosonde pressure and 
temperature measurements tend to be magnirieo by the quadratic terms of (3-2). 
Because the fms radiosonde errors are constant rather than a percentage of 
the measurement, the errors have a greater effect at higher altitudes 
where the measurement values are smaller. For instance, the rms pressure 
error at sea level is about 0.1 percent of the ambient pressure, but at 
15 km altitude the rms error has grown to 1 percent of the ancient 
pressure [21. 

The effects of measurement errors are minimized by using more than 

the minimum of six radiosonde balloons and performing a multiple regression 

to find the refractlvity coefficients N , N„, N,, N„., N,. [3]. In 

r D 9 09 “0 99 

this investigation eight radiosonde profiles were used, providing a 
regression with two degrees of freedom. 

Ifost of the radiosonde data end at a height of about 15 km, and 
refractlvity above the cut-off point must be extrapolated (see Section 4). 
Because of the extrapolation and the greater effects of measurement errors 
at higher altitudes, a four-coefficient model is used at heights above 15 km. 

N = N + 6N. + 9 sin (0) + 99 sin (9) . (3-3) 

r 9 9 99 

This model gives four degrees of freedom in the regres* 'n at the more 
critical higher altitudes with little sacrifice in accuracy, since re- 
fractivlty is more uniform at these heights than it is near the surface. 

Because of balloon malfunctions, only about twenty sets of data with 
eight radiosonde profiles were available. Ten additional sets with seven 
profiles w».re available. However, the quadratic nodel of (3-2) would have 
provided only a single degree of freedom in the regression. To reduce the 
effects of radiosonde measurement errors in these seven-station sets, the 
four-coefficient model (3-3) was used at all altitudes. 



Using the appropriate model, (3-2) or (3-3), multiple regression is 
used to find the coefficients at each height from refractivity measure- 
ments at that height. Once the coefficients have been determined, they 
can be used to calculate the refractivity along a ray in any direction. 

A ray trace which uses a spherically symmetric refractivity profile 
is employed to test the accuracy of the surface correction formula. 

Marini and Ibirray obtained a spherically symmetric profile from measure- 
ments made by a single radiosonde. However, the regression coefficients 
can also be used to obtain a spherically symmetric profile. The profile 
is generated by calculating the refractivity directly above the laser 
site, using the regression coefficients. This approf^h tends to minimize 
the effects of errors in the radiosonde data, and is therefore more 
accurate than simply using the measurements made by a single radiosonde. 
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4. RAY-TRACIIKJ PROCEDURE 

The data used to construct refractlvlty profiles were gathered in 
?r-oj< ct Haven Hop I during January and February of 1970 [41. The Haven 
Hop data consist of measurements of pressure, temperature, and relative 
hunidity made by radiosonde balloons released from sites in the 
Washington, D.C. area (see Figure 2). The balloons were released from 
the 'ites within a few minutes of each other at various times during the 
nigl t and day and tracked to an average altitude of about 15 km. The 
radiosondes report observations of pressure, temperature, and relative 
huffllJity every 30 seconds as they ascend. From these observations, the 
phase refractivLty and group refractivity can be calculated. 

The refractivity must be known at certain standard altitudes between 
tae surface and 1000 km, the assumed satellite height. For altitudes 
b.'.low the radiosonde cut-off point, the radiosonde measurements are 
interpulated t' the nearest standard altitude, following the procedure 
in Ztnr.er, Gardner and Rao [2]. The radiosonde heights are inferred 
from the measurements using the hydrostatic equation. Above the radiosonde 
cut-off, the last measurements must be extrapolated. As is done in 
Zanter *;t al. [21, pressure is assumed to decay exponentially, while 
temperature .«i relative humidity remain constant. Once the pressure, 
temperfrure, and relative humidity are known at the standard altitudes, 
the refractivity can be calculated. 

Some error is introduced because the radiosonde balloons take from 
45 minutes to an hour to ascend, whereas a laser pulse would pass through 
the at:'Ot-,jhere almost instantaneously. This problem can be dealt with in 
e^:ner of two ways: 






( 1 ) 


> 
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It can be argued that the atmosphere does not change significantly 
under normal conditions during ascent, so the ascent time can be 
neglected. In this case, the tracking data are used to determine 
the radiosonde position as it is blovm downrange during ascent 
(tracked ascent). 

(2) Since the balloons ascend on an approximately linear path (that 
is, wind speed and direction are nearly constant at all 
altitudes), the measurements taken along the balloon path are 
a good estimate of the conditions directly above the release 
site at the instant the balloon started its ascent. In this 
case, tracking data are ignored, and the balloons are assumed 
to ascend directly above the tracking site (vertical ascent) . 

Since neither argument is entirely satisfactory, ray traces were made 
under both assumptions. The difference between ray traces on the same 
set of data under the two assumptions was small (see Section 5) . 

At each standard height, a multiple regression is performed, using 

•k 

the eight refractivity measurements and corresponding positions in the 
appropriate refractivity model, (3-2) or (3-3). A three-dimensional 
refractivity profile is thus constructed from the eight radiosonde profiles, 
which allows one to compute the refractivity along any arbitrary ray path. 

A ray path is specified by azimuth and elevation angles from the release 
site. 

Given an azimuth angle and an elevation angle from the release site, 
the regression coefficients and the appropriate model of Section 3 are 
used to generate a refractivity profile along the specified path. This 
profile is then passed to a ray-tracing program, which computes the 
range error. Three types of ray traces are made for each set of data: 

* 

Some sets of data contained only seven radiosonde releases, as explained 
in Section 3. 
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(1) SlnRle Radiosonde Spherically Synnetrlc Ray Trace (RT^) — 

RTjjj uses only the data obtained by the radiosonde released 
from the laser site. Such a ray trace was used by Marini and 
Murray to test the accuracy of their surface correction formula. 
RT^ is made under the assutq>tion of a spherically symmetric 

fBl 

atmosphere. In this report, RT„„ i.<- used to establish agreement 

nn 

with Marini and Murray. 

(2) Spherically symmetric ray trace (RTj^) — Although RT^ is made 
under the assumption of a spherically symmetric atmosphere, the 
use of data from a single radiosonde may cause signlfic-^nt 
errors in the range correction obtained [2]. RT^ uses a 
spherically symmetric refractivlty profile obtained by calcu- 
lating the refractivlty directly above the laser site using the 
three-dimensional regression profile. This approach tends to 
minimize the effects of radiosonde measurement errors (see 
Section 3). Since spherical symmetry is assumed, RTj^ is inde- 
pendent of azimuth, and is not sensitive to the effects of 
horizontal refractivlty gradients. RT^^ is used to test the 
accuracy of the surface correction formula derived in Section . 

(3) Three-dimensional ray trace (RT^) “ RT^ uses the three-dimensional 
regression coefficients to generate a refractivlty profile along 
the ray path. RT^ is dependent on azimuth and contains the 
effects of horizontal refractivlty gradients. Since RT^^ contains 
no gradient effects, the difference RT^ - RTj^ isolates the con- 
tribution of gradients to the range correction. Rays are traced 
every 10° azimuth to find the effects of horizontal refractivlty 
gradients . 




-- i 


Each of the three types of ray traces Is made at four different ele- 
vation angles: 10®, 20°, 40°, and 80°. The centrally located Leonardtovm, 

IM. station (Site 54) w..'. ..sed as the laser site, so that rays could be 
traced at all azimuths. The ray-trace routine employed in this report is 
the Thayer method [5]. This same method was also used by Marini and 
Murray [i] and Zanter, Gardner and Rao [2]. Agreement, with the results 
obtained by Marini and Murray was checked by feeding the radiosonde data 
appearing in Appendix 4 of their report into the ray-tracing program. 

The range correction obtained agreed with that reported by Marini 


and Murray. 


5. RESULTS 


5.1 Horizontal Refr'ictlvlty Gradients 

A spherically symmetric atmosphere would have horizontal gradients 
equal to zero. As stated earlier, nonsymmetric behavior of the atmosphere 
will introduce errors into the surface correction formula and in any ray 
trace made under the assumption of spherical synanetry. 

The horizontal refractlvity gradients can be easily computed from 
the refractlvity model (3-2): 

SN 2 

H = Ng + 4 cos (0) + (0 cos (0) + sin (8)]'i>Ng^ + 20Ngg + T sin (20) 

(5-1) 

|| - sin (0) + 0 sin (0) + 2$ sin^ (0) (5-2) 

where 0 Is proportional to horizontal displacement in a north-south 
direction, and ^ sin 0 is proportional to horizontal displacement in an 
east-west direction. 

Figures 3 and 4 are sample plots of the refractlvity and the north- 
south and east-west gradients versus height. The gradients were calculated 
for points directly above the laser site. The magnitude of the north-south 
gradient is generally larger than that of the east-west gradient. This is 
to be expected, since temperature has a large north-south gradient and 
refractivity is inversely proportional to temperature. For the Haven Hop 
data, surface pressure is approximately constant, while surface temperature 
Increases to the south and east. We would expect refractivity to decrease 
to the south and east, and consequently, both the north-south and east-west 
gradients to be negative. Both figures confirm our expectation. 

In both figures, the north-south gradient reverses sign as height 
increases. The sign reversal is predicted by the hydrostatic equation: 
from Appendix C we may write 









HEIGHT (km) 
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REFRACTIVITY AND GRADIENTS 


Figure 4. Group refractivity N , north-aouth horizontal refractivlty gradient 

3N ® 1 

, and east-west horizontal refractivity gradient 

directly above site 54 at 2330 GMT on 1/21/70. 


H 


20 

/ x-Mg/RB 

' ■ y 

T •= T + $h , 
s 

Approxlaatlng refractivity by the dominant first term of (2-1), we obtain 


P. / 8h -(Mg/R3+1) 
N = 80.343f(X) ^ 1 + f^ 

* s \ S; 


Then taking the derivative in a horizontal direction. 


pt ^1 

N' - N <rr^ + ~ 

g gVs 


-is + lU 1 


(5-3) 


where N', P', T' are horizontal gradients of group refractivity, 
g s s 

surface pressure, and surface temperature. 


If the pressure gradient is ignored, (5-3) predicts a sign reversal at 
h = ^ Tg meters, or about 8 km, which is consistent with Saastamoinen's 
[6] results. 

5.2 Rt /-T racing Procedure 

In Section 4 it was pointed out that the ascent time of the radio- 
sondes introduces an error into the ray-trace calculations, since the 
measurements are not taken at all altitudes simultaneously. One way of 
dealing with this problem is to assume that the meas rements au a point are 
a good estimate of the conditions there at the instant of b'’''loon release 
(tracked ascent). The other is to assume that the measurements are a good 
estimate of conditions directly above the release site, in which case the 
balloon is assumed to ascend vertically. Ray traces for seventeen sets of 
data were made under both assumptions for comparison. The results are 
■presented in Table 1, and a sample comparison of three-dimensional ray 
traces for the same set of data Is shown in Figure 5. There is very little 
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difference between the cwo assumptions. However, since not all balloon 
releases were tracked, some sets of data could only be processed using 
vertical ascent ray traces. Furthermore, Gardner and Hendrickson's 
results appear to agree slightly better with the vertical ascent ray 
traces [71. Therefore, the remaining results presented used vertical 
ascent ray traces. 


TABLE 1. 

COMPARISON OF SPHERICALLY SYMMETRIC RAY TRACES 
USING tracked ascent 'UiD VERTIC.\L ASCENT (17 SETS) 


Elevation 

RTj^ (vertical) - RTj^ (tracked) 
Mean (cm) SD (cm) 

10“ 

.29 

.75 

20“ 

.11 

.46 

40“ 

.06 

.25 

30“ 

.04 

.16 


The surface correction formula derived by Marini ana Murray [1] (see 

Appendix D) was designed to agree with the single radiosonde spherically 

symmetric ray trace RT,_,. However, as seen in Table 2, a bias exists 

between RTj^ and Since RTj^ uses the regression coefficients and is 

thus less sensitive to radiosonde measurement errors than RT^„, the surface 

nn 

co'.rection formula developed in Section 2 was designed to agree with RT^. 

The error in Marini and Murray's fo'*mula (MM) is given bv MM - RT.„,, while 

Mn 

the error in the formula of Section 2 (AR) is given by - RT^^. These 
errors are compared in Table 2. While the errors in both formulas have 
about the same itandard deviat i.in , t'r. formula of Section 2 has a near 
zero mean, whil Marini and Murray's f irmula • 'uS. This bias 


was also observed by Zanter, Gardner ai 




TABLE 2 


C(»fPARlSON OF SPHERICALLY SYMIETRIC RAY TRACES 



^^1 “ ''''mm 

(24 sets) 

MM-RT^ 
(24 sets) 

AR - RTj^ 
(31 sets) 

Elevation 

Mean (cm) 

SD (cm) 

Mean (cm) 

SD (cm) 

Mean (cm) 

SD (cm) 

10* 

0.54 

0.75 

-0.40 

0.49 

-0.03 

0.46 

20“ 

0.29 

0.40 

-0.25 

0.25 

0.06 

0.25 

o 

0 



0.15 

0.21 

-0.20 

0.14 

0.00 

0.14 

80“ 

0.10 

0.13 

-0.13 

0.09 

0.01 

0.09 


RT^ - spherically symecrlc ray trace range correction 
RT^ > single radiosonde spherically synmetric ray trace range 
correction 

IM * range correction predicted by Marini and Murray's formula 
AR * range correction predicted by surface correction forimila 
in Section 2. 

Ten of the thirty-one available sets of data contained only seven 
radiosonde releases rather than eight, and these sets were processed using 
the four-coefficient refractivity model [Eq. (3-3)] to reduce the error in 
the regression coefficients. To examine the differences between ray traces 
processed using the six-coefficient model [Eq. (3-2) J and those using the 
four- coefficient model, another ten sets of data containing eight radio- 
sonde releases were processed twice — once using the six-coefficient model 
and once using the four-coefficient model. The results are summarized in 
Table 3, and a sample plot of RT^ vs. azimuth for the same set of data using 
each model appears in Figure 6. 

The four-coefficient model shows very close agreement with the six- 
coefficient model, particularly in the means, which are near zero. The 
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AZIMUTH 

Figure 6. Comparison of RT- for ray traces using six-coerf iclent model (Eq. 3-2) and ^ 

four coefficient model (Eq. 3-3) at 10“ elevation. Balloons released at 
1530 GMT on 2/16/70. 
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four-coefficient ray traces agree extremely closely in gradient, effects 
(RT^ - RTj^) with the six-coefficient model. The extra quadratic terms of 
Eq. (3-2) apparently contain little information about the horizontal 
gradients, and so the predominant gradient effects appear to be linear. 


TABLE 3. 

CO{PARIS(»l OF SIX-COEFFICIENT AND FOUR-COEFFICIENT MODELS (10 SETS) 


RTj(6) - RTj(4) 


RT3(6) - RT3(4) 


(RT3 - RT3H6) - 
[RTj - RT^](4) 


Elevation 

10 " 

2C* 

40“ 

80" 

No. of 
Obs . 


Mean (cm) 


0.09 


0.05 


0.03 


0.02 


SD (cm) 


0.85 


0.43 


0.23 


0.15 


Mean (cm) 


0.00 


0.04 


0 . 0 ’ 


0.02 


SD (cm) 


0.79 


0.40 


0.21 


0.14 


Mean (cm) 


-0.09 


- 0.01 


0.00 


0.00 


SD (cm) 


0.29 


0.06 


0.01 


0.01 


10 


360 


360 


5.3 Three-Dimensional Ray Traces 

For each set of data, a three-dimensional ray trace RT^ Is made at 
every 10 “ azimuth from the laser site for each of four elevation angles: 

E ■ 10“ , 20“ , 40“ , 80“ . The effects of horizontal refractlvity gradients 
are isolated by subtracting the spherically symmetric range error RT^^, 
which is independent of azimuth, from the three-dimensional range error 
RT 3 at each azimuth angle. Sample plots of RT^ - RT^ vs. azimuth appear 
in Figures 7 and 8 . 

The E = 80 “ curve shows the effects of computer round-off noise, since 
the refractlvity differences computed in the ray trace routine are quite 
small. Except at E - 80 “ , the range error is an approximately sinusoidal 
function of azimuth, having a maximum towards the south and a minimum 
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0® 30® 60® 90® 120® eo® ISO® 210® 240® 2TO® 300® 330® 360® 

AZIMUTH 

Figure 7. Difference between three-dimensionar ray trace (RT_) and spherically syrametric 

ray trace (RT,) at lO” and 20® elevation’. Balloons released t 1130 GMT on 2/16/70. 
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Figure 8 


Difference between three-dimensional ray trace (RT«) and spherically symmetric 
ray trace (RT^^) at AO® and 80® elevation. Balloons released at 1130 GMT on 2/16/70. 
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towards the north. One might expect just the opposite, since at the surface 
there are colder temperatures, and therefore larger values of refractivity 
to the north. If the north-south refractivity gradient remained negative 
at all altitudes, one would expect to trace through denser air to the north 
than to the south, and consequently to have a larger range error to the 
north. However, as noted earlier, the north-south gradient changes sign 
around 8 km. The positive gradients above 8 km overcome the effects of 
the negative gradients near the surface, because the refractivity 
differences are smaller near the release site than they are at higher 
altitudes further from the site. Thus, even though the refractivity 
itself is large at lower altitudes, the contribution to the range error is 
from refractivity differences, which are small. 

The means and standard deviations of RT^ - RTj^ combined over all 
azimuths for the thirty-one sets of data appear in Table 4. The means 
are all very close to zero, and so the difference RT^ - RT^^ is a good 
measure of the effects of horizontal refractivity gradients. The effects 
of these gradients are quite significanL at the lower elevation angles. 

The means and standard deviations of RT^ - RT^^ for each azimuth angle are 
plotted in Figures 9-12. At E = 10®, the mean swings ± 2.5 cm, but the 
swing drops to less than ± 0.2 mm at E = 80°. The mean curves are very 
close to being sinusoidal, except at E = 80°. The standard deviation 
appears to be a function of azimuth. However, this effect is due to the 
limited number of surface stations available (seven or eight) and the 
regression used in performing the ray trace. The regression error is 
discussed in detail in Gardner and Hendrickson's report [7]. 
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Figure 12. Standard deviation of difference RT» - RT. versus azimuth at 
40® and 80® elevation (31 sets of data). 





TABLE 4 


RESULTS OF THREE-DIMENSIONAL RAY TRACES 
(31 sets, 1116 observations) 



RTj - RT^ 

Elev. Angle 

Mean 

IE 

10° 

-0.0666 cm 

1.9009 cm 

20° 

-0.0079 

0.4968 

O 

O 

-0.0008 

0.1160 

80° 

-0.0005 

0.0127 


Histograms of RT^ - RT^^ for the four elevation angles are plotted in 
Figures 13-16. Since RT^ - RTj^ is approximately sinusoidal, one would 
expect the histograms to resemble the aenslty function for a sinusoid with 
random phase and amplitude. For a fixed amplitude A and random phase 6 
uniformly distributed between 0 and 2ir, the probability density Y = A sin 6 
is given by 


r 


f.Cy) = ( 


,Va 2 - y2 


Vv 


y| < A 


otherwise 


f^ has peaks at ± A and a minimum at y = 0. However, the random amplitude A 
tends to smear the peaks. The smearing nearly obscures the peaks in the 
E * lO" histogram, but they are apparent at E = 20° and E = 40°, where the 
amplitude varies less. The E = 80° histogram contains computer round-off 


noise. 
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Figure 14. Histogram of difference RT^ - RT at 20® elevation 
(31 sets of data, 1116 observations). 
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6. CONCLUSION 

The effects of horizontal refractivity gradients on the accuracy of 
laser ranging systems was investigated by ray tracing through three-dimensional 
refractivity profiles. The profiles were generated by performing a multiple 
regression on measurements from seven or eight radiosondes. The refractivity 
models provided for both linear and quadratic variations in the horizontal 
direction. To isolate the gradient effects the range calculated from a 
spherically symmetric ray trace (RTj^) was subtracted from the range calculated 
from the three-dimensional ray traces (RT^). The mean of the difference 
RT^ - RTj^ was a sinusoidal function of azimuth having a minimum near 0* 
azimuth (due north) and a maximum near 180** azimuth (due south). The peak- 
uo-peak variation was approximately 5 centimeters at 10“ elevation and 
decreased to about 3 millimeters at 40“ elevation. The standard deviation 
of RT^ - RTj^ was approximately 1 centimeter at 10“ elevation and decreased 
to 0.5 millimeters at 40“ elevation. The ray trace results also indicated 
that the linear variation of the refractivity in the horizontal direction is 
the primary error source in the range correction formulas. The effects of 
quadratic and higher-order variations appear to be negligible. 
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APPENDIX A. 


APPROXIMATION OF 


sin (6) 


sin (6) sin (0 q) 




tan (6 q) 


Consider 1 - 




- i voj2 . 1 ^ lo-v^r — y 
\ L(1 + 10 ^N)(r» + h). 


= 1 - ll + 10'^N_V 7 r \ T — ^ 

' °/ 1 + fl0~^N + — + 10"^N — 


' 1 - + 2 X lO’^Ny + 10 "^^NqJ 

• I (i + l)f-io'S - ^ . 

i=0 ^ 0 0 ' 


,-12.,2 


(A-1) 


Neglecting 10 Nq (Nq = 300) and taking the first two terms of the series 




= 1 - 1 + 2 X 10 


\)(' - 


2 X io~^N - — - 2 X 10"^N — 

’’o *^0 


. 


Taking only first-order terms. 


1 - (-^1 ^ 2 — - 10 ^(N - N) 

\ ^0 ° J 


, X ^ 3X ^ ^ , 

1 - -;r + for X - ' I 
2 8 


Equation (A-1) can be written 


liJi 







APPENDIX B. 


SIMPLIFICATION OF EQUATION (2-9) 


Consider the —— terms of Eq. (2-9) 


10 (Nq - N) Ng dh 


+ I • 10“^ 


10 " 
sin (E) 


10’^^ f 


I ai> + io-‘ f - io-Nn„ - N 

■ ^ 

fl - 10'*S^ + I • I0‘*^(llj - N)^ I N 

E) fl - 3 " I "“g ’*'• 

3 10-6 

f 2 • 10 f 2 

1 NN dh + h^N dh 

J 8 ,-2 rpi J 8 


10 (N. - N) N dh 

0 g 


+ ^7 — 7^ 1 - 3 X 10"^(N 

r^j sin (E) ' i 


hN dh 
g 


+ -7 — 77T 1 NN dh + 
sin (E) J g 2 

* A 


I . 10-^ 


sin (E) 


in~6 f in~6 f in~^2 

1 - 10 N- N dh + 7 ^, hN dh + NN dh 

sin (E) 0 J g r_ sin (E) J g sin (E) g 

1 . in~6 

2 10 f 2 

+ h N dh 

r- sin (E) ■' ® 


where N^ = 300 ^ N, E 10“, and the evaluations [Eq. (2-11)] have been 
used to eliminate any terms contributing less than 1 mm. A similar 


analysis of the remaining terras gives 


AR = (lO ^ f N dh + lO"^ I (— - 10“^(N„ - N)1 N dl)» 

sin (E) j g J 0 j g j 


-6 r f'h 

i ’’n 


i dl^ 

^ ; 


- 4 — ^ hN dh 

sin^ (E) I J ® 


+ 10 [ (NN - J N^)dh + — ! h^N dh) 

J g 2 2 J g j 


1 f 3 10"^ [ ,2., .. . 3 ' lO"'^ ! 

+ p \ 7 * — :r~ h N dh + ; hNN dh 

sln^ (E) V J « 'O ^ . 


Consider the following pair of terms of (B-1) 



3 X 10~^ 
Tq sin^ (E) 


h^N dh 
g 


-12 

3 X 10 
Tq sin^ (E) 


I 


hNN dh 
g 


-1.64 X 1 0 
• 3 

sin (E) 

6.92 X io~ 
sin^ (E) 


42 


where the evaluations of (2-11) or Appendix C have been used, with 
Nq - 300, Pg = 1000 mb, = 273"* K, and K == 0.9. 

At E = 10®, the sum of these terms is about 1 mm, and they approximately 
cancel. At E ^ 20®, each of the te:ms is individually less than 1 mm and 
they are thus insignificant. For 10® < E < 20® the sum of these terms is 
1.00 mm or less. Leaving these terms out, (2-10) is obtained. 
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APPENDIX C. 
INTEGRAL EVALUATIONS 


1 . 


h^N dh 
g 


Approximating N by the dominant first term of (2-1), 
g 


h N dh = 
g 


h^(80.343)f(X) I dh . 


Marini and Murray use the hydrostatic equation to obtain [1, Appendix 2] 


P = P 


_ 

(T_\ 6R 




and T = T + 6h 
s 


where M = 28.966 = moleculfr weight of dry air (kg) 


g = 9.8 = acceleration of gravity (m/s) 


j oules \ 


R = 8314.36 = universal gas constant - kg mol ej 
8 = temperature lapse rate 


P^,T^ = surface pressure and temperature. 


h^Ng dh == 80.343f(A) ~ 


r T + 6h 


3R 


+1 


h dh 


s / 


Integrating by parts twice, we obtain 


, . 3 2 

9 / R \ 2 9K 

h N dh 80.343f(A) P T ^ - 

g \Mgj s s 2 - K 


where K = 


1 _ M 


Mg 


and so 


10 


-6 


h N dh 

g 


1.4961 
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2 . 


hNN dh 
g 


Let N = Ng = 80.343f(X) ^ . 


hNN dh = 
g 


(80.343)^f(X) ^ h dh 


Integrating by parts, 


hNN dh = (80.343)^f(A) ^ , 

g BT ^ ^ J 

R0 


r n + 6h\ 




dh 


= (80.343)^f(A)P^I-r^l ^ — - 


s^2Mgl 3 - 1/K • 


Thus we obtain 


3 X 10 


-12 


0 


hNN dh = 6.533 x 10"^^P^ ^ 


s 3 - 1/K 


10 


-12 


(N_ - N)N dh 
0 g 


This integral divides into two terms which have been evaluated in [1] 
■12 


10 


(N^ - N)N dh 
0 g 


= 10 ^Nq^IO"^ Ng dh^ - 10 


NN dh 
g 


10 


.002357P 


f(A) 
s f(4>,H) 


2 X 4.7343 X 10"®f(A) '' ^ 


T 3 - 1/K 
s 


where the less significant terms have been ignored. Approximating N : 

J 


~ 80.343f(A) ^ 
s 


we obtain 





lO'^^ I (N_ - N)N dh = 9.4682 x io~® f(X) 
J 0 g T 


i 


' •‘*%tSWK®rrJBS,iW‘ 
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APPENDIX D. 

MARINI AND MURRAY'S SURFACE CORRECTION FOR>IULA 


Marini and Murray [1] developed a surface correction formula using an 
analysis like taat in Section 2. They considered only the three most 
significant terms of Equation 2-10 (the underlined terms), and obtained 


MM = 


f (A) A -f B 

■ sin (E) + 


sin (E) + .01 


where A = 0.002357P + 0.000141e 

s s 

B = 1.0R42 X 10’®P T K + 4.7343 x 10~® -- - — 

SS Xj'"X/K 

S 


and f(A), F(6,H) and K are as defined in Section 2. 

The starred B terms are an "optional adjustment" to reduce a bias at 
elevations near 90°. However, our investigation showed these terms to have 
a significant effect at all elevations. Including the starred B terms 
increases the value of MM by about 1.6 cm at 10° elevation, 0.8 cm at 
20° elevation, 0.45 cm at 40° elevation, and 0.3 cm at 80° elevation. 

The starred B terms were neglected in our comparisons using Marini and 


Murray's formula. 
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